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The path integral, which generates in-in correlation functions of a scalar field in a cosmological 
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of the instantons in Euclidean field theories. When the scalar potential has more than one locally 
stable vacua, the correlation functions receive contributions from all of them via these instanton-like 
configurations, which is similar to tunneling. We present some explicit solutions for toy models and 
discuss possible implications of our results. 
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I. INTRODUCTION 

Understanding the physics of the early universe is both a challenge and an opportunity. With the recent obser- 
vational advances in cosmology, one may hope to test new theories on energy scales that can never be reached by 
conventional accelerators and discover signs of various novel ideas. This also motivates one to have a comprehensive 
understanding of conventional approaches. 

Inflation is the most attractive paradigm in solving the problems of the standard big-bang model and its main 
predictions are in agreement with observations so far. However, most of these predictions are based on semiclassical 
reasoning and a better understanding of inflationary theories in terms of well established physics is necessary for 
their robustness. In recent years, following the work [1, 2], there has been a growing interest in calculating quantum 
loop corrections to cosmological correlation functions during inflation. Earlier, various cosmological implications of 
(loop) quantum effects have been studied (sec e.g. [3-15]). One of the main reasons for the recent interest is the 
observation that interactions give rise to primordial non-gaussianities [16], which can potentially be observed (see [17] 
for a review). 

As usual, the main technique in evaluating quantum corrections is the (in-in) perturbation theory. Perturbative 
loop calculations corresponding to massless fields in inflationary spacetimes are plagued by infrared divergences and 
there are different views in the literature whether infrared effects are real or not (see e.g. [18-32], see also [33-35] for 
earlier work) . On the other hand, loop contributions may also be anomalously sensitive to the UV cutoff as discussed 
in [36]. There are some attempts to get non-perturbative results such as the stochastic approach [37-41]. There are 
also some attempts to go beyond the one loop approximation [42-44] or to work out the complete one-loop effective 
action in a time dependent background [45-51]. 

In this paper, we consider the in-in path integral for the generating functional of cosmological correlations corre- 
sponding to an interacting scalar field in a Freedman-Robertson- Walker (FRW) spacetime. For a field theory defined 
in the flat space, instantons, if they are admitted, give valuable non-perturbative information about the vacuum 
structure of the theory (or one may say that instantons exist if vacuum is nontrivial). Our aim here is to see whether 
similar instanton-like classical solutions exist, which would correspond to nontrivial saddle points of the cosmological 
in-in path integrals. In searching for instanton solutions, one naturally makes a Wick rotation to the Euclidean sig- 
nature [52]. However, it is not known how to Wick rotate a general cosmological spacetime. As an exception, if one 
considers the de Sitter space in conformal coordinates, which is appropriate for Wick rotation due to the special form 
of the metric, one sees that Euclidean continuation does not give an "inverted" scalar potential. Therefore, efforts of 
constructing Euclidean solutions directly reminiscent of the usual instantons would fail. 1 

On the other hand, a cosmological in-in path integral is very different in a few important ways from an in-out path 
integral defined in flat space. Namely; 

(1) The boundary conditions arc different compared to the in-out path integration. 

(2) In applying saddle point approximation to an in-out path integral, one searches for classical configurations which 
have finite action to have a well defined expansion scheme. However, in-in path integral contains two action terms 
in the exponential with different signs and there is a possibility of cancelation even for configurations with "infinite" 
action. 

(3) The cosmic expansion changes the classical dynamics in a crucial way. 

Thus, one may still search for classical configurations corresponding to stationary phases of an in-in path integral in 
the Lorentzian signature. As we will discuss, all the three properties listed above will be important in applying a 
rigorous stationary phase approximation. 

Indeed, we will be able to show that when the scalar potential have more than one locally stable vacua, there are 
nontrivial classical configurations which correspond to stationary phases of the in-in path integral for the generating 
functional. Due to the existence of these solutions, the correlation function receives contributions from all of the 
vacua. For instance, the vacuum expectation value of the scalar (1-point function) becomes the sum of the all of the 
minima of the potential even when the in- vacuum is specified around one of the vacua, which is like a tunneling effect. 

As we will discuss, when the in-vacuum state is defined not at the asymptotic past infinity but at a finite time in the 
past, there appears some additional technical difficulties mainly related to the identification of the vacuum and the 
correct set of boundary conditions one must employ. We show that with some reasonable assumptions it is possible 
to construct instanton-like solutions for that case also. 

The plan of the paper is as follows. In the next section, we review the in-in path integral formalism for an interacting 
scalar field propagating in a cosmological spacetime. Most of the material presented in that section is well known. 
As partially new results, we give a slightly different derivation of the equivalence of the path integral and operator 



1 Here, we only consider quantum fields in a fixed background. There exits gravitational instantons in the context of Euclidean quantum 
gravity in three dimensional de Sitter space [53] . 
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formalisms and we discuss Wick rotation for a scalar in de Sitter space (see also [54] for a rigorous discussion of 
Wick rotation). In section III, we fix boundary conditions and apply stationary phase approximation to determine 
the equations obeyed by saddle points. It turns out that the results alter if the in-vacuum state is define at finite or 
infinite times. In that section, we give several examples and elaborate physically on why instanton-like solutions do 
not exist in some situations and why they are admitted in some others. In conclusions, we discuss possible implications 
of our results and indicate some open problems. 



II. REVIEW OF THE IN-IN PATH INTEGRAL FORMALISM 



In this section, we would like to review the path integral derivation of the in-in correlation functions. Our aim is 
to fix our notation and discuss Wick rotation to the Euclidean signature as a first attempt to search for instanton 
solutions in dc Sitter space. The background metric is taken as 

ds 2 = -dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ) 

= a(r,) 2 (-di 1 2 + dx 2 +dy 2 + dz 2 ), (1) 

where t and r\ are the proper and the conformal time coordinates, respectively. The scale factor of the de Sitter 
space is given by a = cxp(Ht) or a = —i]o/r], where rjo = 1/H. Most of our considerations below will be valid for an 
arbitrary scale factor a(t) and we will freely switch between these two coordinate systems. The action of a minimally 
coupled real scalar field propagating in this background is given by 

S[^J] = - l -j d^^S [V„0 V"0 + V (<£)] - 2J<j>, (2) 

where we introduce an external source J, which can be turned on or off. In this paper we only consider real scalar 
fields with this canonical action. 



A. General theory 

Our main object of interest is to calculate the vacuum expectation value 

<0,t |Q(i)|0,io>, (3) 

where Q(t) is a polynomial of the field variables and |0, t > is the ground state of the system at time t - Preferably, one 
would let to — > — oo, however in a realistic scenario to can be finite as well (in that case there is a certain uncertainty 
even in the specification of the free vacuum [55]). One should also work with the vacuum of the interacting theory 
and thus in perturbation theory corrections to the initially chosen free vacuum state must be taken into account. 
The most straightforward way to introduce a path integral for (3) is to consider the generating functional 

Z[J+,J-}= f D<l><0,to\<t>,t>j-«t>,t\^to>j+, (4) 

where / D<j>\<j), t><<p, t\ is the identity operator constructed from the field variables at time t and the inner products 
in (4) are evaluated in the presence of two independent external sources J + and J~ coupled to field variabes. Differ- 
entiating (4) with respect to J + or J~ at time t, and setting J + = J~ = gives (3). It is easy to write the transition 
amplitudes in (4) in terms of path integrals to obtain 

Z[J+,J-} = j ' D4>j ni?0+2?re lS[0+ ' J+1 - j5[ ^'^ ] *o[0 + (io)]*5[r(io)], (5) 

to 

where <j) + and <j)~ integrals are over all field configurations starting from time to and ending at time t obeying 

(j) + {t,x) = (l)-{t,x)=cj>{x), (6) 

and \E , o[0 ± (*o)] are the vacuum wave-functionals corresponding to the inner products <0 ± (t o )|O, t >. By introducing 
a Dirac-delta functional it is also possible to rewrite the path integral as 

Z[J+, J-] = f n^+2?0- e jS ^ + ' J+ ]- jS [^'^]* o [0+(t o )]*S[^(to)]<5[0 + (t) - r(t)}, (7) 

J t 

where now there is no restriction imposed, and (7) reduces to (5) after integrating over 4> + (t) or (p~(t). 
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B. Perturbative expansion 



Before discussing the Wick-rotation in Sitter space, it is instructive to see how perturbation theory works in the 
path integral formalism. 2 For that the free scalar action can be written as 



So[4>, J } = \f d i x[cj>L ( j>]+2Jcj>, 



(8) 



where L is a (second order) differential operator. It is known that the basic role of the vacuum wave-functionals in a 
(free) path integral is to impose the proper ie prescription necessary to define a unique inverse of L (see section 9.2 
of [56]). Therefore, using an integral representation of the delta functional in (7), the free generating functional can 
be written as 

Z [J+,J-} = J DX(x) J JjD^ + D0-e iSo ^ + - 7+ ]- <s °[*" >J "]exp (i J d 3 xX(x) [(j> + (x,t) - <p~(x,t)] 

to 

= J DX(x) J J|2?$expQ^ dt'd 3 x [$ T L$ + $ T j] 



(9) 



where we define 3 



L 
—L 



J 



J+(x,t') + X(x)S(t' -t) 
-J-(x, t') - X{x)5{t' - t) 



Performing the Gaussian integral over $ one finds 



Z [J + ,J-}= J DX{x)cxp(^-^ £ dt'd 3 x' dt"d 3 x"J T AJ^j , 



(10) 



(11) 



where A is the inverse of L (with the ie prescription implied by the vacuum wave-functionals 4 ). One can write A as 



A++ A+- 
A-+ A— 



(12) 



where LA ++ = 1, LA = — 1, LA+ = and LA + = 0. Since L is a symmetric operator Green functions obey 

A+- (f , t'; x", t") = A~ + (x", t"; x' , t'), 

A++(x',t';x",t") = A + +{x\t"-x',t t ), (13) 
A (x r , t'- x", t") = A (x", t"; x',t'). 

Inside the exponential of the functional integral (11), there are quadratic and linear A terms. After performing the 
integrals over t' and t", the quadratic terms become 



exp ^-i J d?x'd 3 x"X(x') [A++(t) + A (t) - A+~{t) - A~+{t)] A(x")) • 



(14) 



The Green function inside the square bracket is annihilated by L and thus it is not invertible. Therefore, the only 
way to make A path integral to be well defined is to impose these terms to cancel each other. In that case A integral 
becomes 

J DAcxp ^-i J d 3 x / d 3 x' / X(x') [[A ++ (i) - A-+(t)]J+ + [A (t) - A+-(t)]J-]^ , (15) 

which is non-zero for arbitrary J + and J~ provided 

A++(i) = A-+(t), 

A"-(i) = A+-(i). (16) 



2 This issue is discussed in the appendix of [1], Our treatment of the delta functional in (7) is different than [1] but its effect turns out 
to be the same, namely to force the appropriate boundary conditions on the propagators consistent with the operator formalism. 

3 Note that it is possible to shift the upper limit of the time integral in (9) by an arbitrary positive number since in-in formalism guarantees 
that the contributions for times larger than ( cancel. This justifies the delta functions in (10). We would like to thank the anonymous 
referee for pointing out this to us. 

4 Although in the free theory this procedure is relatively easy to implement, applying ie prescription for the interacting case using the 
in-in perturbation theory is highly non-trivial. Sec [57, 58] for clarification of this point. 
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These are the boundary conditions, which must be imposed on the homogenous solutions A + and A + to get a well 
defined path integral. When these are imposed A integral decouples and the generating functional becomes 

Z Q [J+, J-} = cxp dt ' d3x ' J t dt"d 3 x" J T Ajj , (17) 

where 

(18) 

Differentiating this final expression with respect to J + and J~ one may verify from the definition of Z that 

A++(f ,t';x",t") = -i <0,t \T<j>(x',t')(l)(x",t")\0,to>, 

A—(x',t';x",t") = -i <0,t \Tcj){x',t')<t>(x",t")\0,to>, (19) 
A" + (f ,t';x",t") = -i <0,t \4>{x ,t')<f>(x" 7 t")\0,t o >, 

where T and T deontes time and anti-time orderings, respectively. These are exactly the propagators obtained in the 
operator formalism. As usual, for an interacting theory with a polynomial potential V((j)), the generating functional 
can be expressed as 

Z[J + ,.r] = exp 
which can be evaluated order by order using perturbation theory. 




-iV(-iS/SJ+)/2 



exp 



i / V(iS/8J-)/2 



z [j+,j-}, 



(20) 



C. de Sitter space and Wick rotation 



At this moment, it is useful to discuss the specification of the vacuum state |0, to>, which is closely related to Wick 
rotation. For that discussion, we specifically consider a real scalar field in the Poincare patch of the de Sitter space and 
use conformal coordinates (rj,x), where r\ < 0. The easiest way to specify the interacting vacuum at 770 = —00 is to 
employ a projection by giving a small imaginary part to the time parameter 77. In general, to project out an arbitrary 
ket- vector onto the interacting vacuum state defined at 770 = —00, one may introduce the operator exp(— eHAi]) 
with A?7 — > 00 and e > 0, where H is the exact Hamiltonian. 5 Since the unitary time evolution operator is given 
by U — cxp(-iHrj), such a projection can be naturally incorporated by assuming that time parameter has a small 
negative imaginary piece. Therefore for fields evolving forward in time, which correspond to the + branch in the in-in 
formalism, the time must be complexified as 



77+ = jy r (l - ie) 



(21) 



where ?7 r is real and e > 0. Similarly, to project out an arbitrary bra- vector onto the interacting vacuum, time must 
be complexified as 



77- = rj r (l + ie). 



(22) 



After these deformations of the integration contours, there is no need to keep the vacuum wave-functionals in (5) and 
one writes 



Z[J+,J-} = f D<t> f []^ + / Y[V<t>-e lS[ * + > J+ 



]-iS[<j,~ ,J~] 



where C+ and C_ are defined by (21) and (22) in the complex Ty-plane (see Fig. 1). 

Let us determine how introducing this complex tilt affects the propagators of a massless field in de Sitter space. It 
is well known that the scalar field operator can be expanded in terms of the mode functions as 



d 3 k 



(2tt) 3 / 2 V2k 



^ik.x — ikrj 






kr] 



x-\-ikr) 
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^ kr] 





(23) 



In the following discussion we take H to be time independent. With slight modifications, the arguments must be generalized to 
time-dependent situations, at least in the adiabatic limit. 



C- 



FIG. 1. The integration contours C+ and C_ in the complex ?)-plane. The Wick rotated contours are also shown as C E 



In evaluating the field operator <f>(x, rj) on the contours C + and C_ , one should replace 77 by i] + and , respectively, 
which gives 



A-+( V '_;v+) = J d 3 ke-^'-+ ik <... = J d 3 ke k <^+*"\.. = J 



d 3 ke~ ke . 



(24) 



and 



& ++ (v+,V+) =0W-V") J d 3 ke~ ikr] + +ikr] '+...+ 6/(77" -rf) J d 3 k e ~ ik < +ik < ... 

= %W) / d 3 ke k «-<\.. + e(r 1 " J d 3 ke k «~^\.. = J d 3 k t 



-kc 



(25) 



where the signs are fixed by paying attention to the orderings along the respective contours and any positive number 
multiplying e is absorbed in e. The crucial point is that properly complexifying the time coordinates gives the extra 
exp(— ke) factors in the momentum integrals, which are necessary for the UV convergence of the Green functions (see 
e.g. [7]). A similar calculation shows that the same damping factors appear for and A , which become well 

defined in the UV. 

The above discussion clearly indicates how one should perform Wick rotation to the Euclidean signature. 6 Since 
the propagators are well defined along the contours C + and C_, and they are analytic functions of their complex 
variables provided e > 0, one can do the following replacements without changing the propagators: 7 

V+ -> -iVr, 
?7- -> +ir) r . 

After this rotation, the Lorentzian actions for + and — branches 



will be transformed into 



zJc ± 



dijd 3 x 



±\2 



dr\ r d x 



(d Vr 



Vr 



2 ' 



(26) 



(27) 



Note that the relative signs of the kinetic and the potential energy terms are the same for the Lorentzian and the 
Euclidean actions since both signs are changed. Of course, this could be anticipated from the beginning since the sign 



In flat space, Wick rotation of the in-in path integrals has been discussed in [59]. 

In other words, by Wick rotation the oscillating functions damped by a small convergence factor are replaced by the exponentially 
decaying functions of the Euclidean time. 
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change for the kinetic term arises due to time derivatives and the sign of the potential is changed due to the metric 
fuctions. However, it was necessary to work out the intermediate steps as we did above to make sure that analytical 
continuation can be performed without a problem. 

As a result, we see that Wick rotation to the Euclidean signature in the Poincare patch of the de Sitter space does 
not give an inverted potential and therefore it is not possible to construct standard instanton solutions for the scalar 
fields, which are supposed to extrapolate in between different vacua as in the case of flat spacetime. 8 Nevertheless, 
we will see in the next section that it is possible to construct classical solutions in the Lorentzian signature, which 
are reminiscent of instantons. In the de Sitter space, if the scalar field is assumed to depend only on time, as we 
will suppose in the next section, then the Lorentzian and the Euclidean field equations will be the same due to the 
special property of the Wick rotation discussed above. Namely any Lorentzian solution would also satisfy Euclidean 
field equations and vice versa. Therefore, in de Sitter space these solutions also correspond to the saddle points of 
Euclidean path integrals similar to instantons. 

It is useful to remind that our discussion in this subsection is carried out in the Poincare patch of the de Sitter space, 
which is relevant for cosmology. The full de Sitter space has different properties and viewing it as the hyperboloid in 
the Minkowski space of one higher dimension, the Wick rotation to the Euclidean signature gives a sphere. 



III. STATIONARY PHASES AND INSTANTON-LIKE STATES 



Consider the following path integral 

/ = J D4>J J]2ty+2ty- cxp (iS[<f>+] - iS[<f>~]) M<P + (to)M[r (to)} , (28) 

*0 

which appears in the generating functional (5). Our aim is to see if stationary phase approximation can be used to 
evaluate (28). We take the scalar action (2) defined in a general FRW spacetime with the metric (1). There are three 
integration variables in (28) and the phase must be stationary with respect to each of them. Moreover, although their 
existence only affects the integral at to, the presence of the vacuum wave-functionals must also be taken into account. 



A. Boundary conditions 

Surface terms may arise in our discussion for two different reasons: cither as a result of integrating by parts the 
field variables or from the variations of the action. The vanishing of the surface terms for the first case is important 
to have a well defined path integral. For example, in the free theory these integration by parts are necessary to obtain 
the second order differential operator in the action. Boundary conditions making these surface terms to vanish must 
be imposed on all paths contributing to the in-in path integral, defining the function space in which the functional 
integral is carried on. The surface terms in the second case will be important in searching for stationary phases of 
the path integral as we will discuss in the next subsection. 

Although it might be important to determine the precise boundary conditions in a more detailed study, here we 
simply assume that necessary conditions are imposed at spatial infinity and concentrate on the surface terms arising 
in the time direction. These surface terms can only arise from the kinetic term in (2). Here, one has the option of 
choosing two different alternatives. If one insists on freely integrating by parts the + and — branches, the following 
conditions must be imposed on the fields separately 



Strong conditions: -g- (t) = 0, ^(t ) ^0)^7^0) = 0. (29) 



On the other hand, in many cases it would be enough to require the absence of surface terms for simultaneous 
integration by parts of the + and — fields, which implies 

6+ 

dt> U;y " v " u/ dt> 



Weak conditions: ^^T^ ~ "g^T^ = °' a ( f o) 



+ (*o&(to)-^-(to)^r(to) 



= 0. (30) 



As noted above, these boundary conditions determine the function space on which the path integration is carried out. 
Our main results will not depend on the choice of the strong or the weak boundary conditions. However, it would be 
important to find out the correct conditions for some applications. 



Note, however, that the relative signs of the kinetic and the spatial gradient terms are changed. 
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To continue, one must consider the infinite and the finite to cases separately. For to finite, we assume that a(to) 
is well defined, i.e. we exclude the situations involving a big-bang singularity. For the other case t = — °°: we will 
assume, having of course inflation in mind, that a(to) — > as to ~~ ^ — 00 (which we simply write as a(— 00) = 0). Let 
us first discuss to = —00 case for which the boundary conditions can be determined unambiguously. We would like 
to remind the reader that the letter t is reserved to denote the present time and we use t' as a dummy variable if 
necessary. 



B. The case t — 00 



The path integral (28) is over all paths extending from t = —00 to time t (+ branch) and then back to to = —00 
(— branch). Therefore, a stationary phase of (28) is a path, conveniently named as $ c ;(t',x), which has independent 
+ and — branches, denoted by cf> cl (t',x), respectively (the branches are connected at time t). The variation around 
such a path can be parametrized by 5<j) + and 5<j)~ obeying 



54>+{t,x) = 6(/>-(t,x) = 54>{x), 
where 8<t>{x) corresponds to the variation of the path at the "boundary" time t. The path $ c ; is stationary if 



(31) 



jZz - S[4T 



0. 



(32) 



From this variation, the following surface term arises (as usual surface terms along spatial directions are assumed to 
vanish by suitable boundary conditions) 



lim 

to — >— 00 



6(j) + a 



dt' 



6(f)- a 3 



dt' 



8<t> 



to 



dt 



+ 

fit) 



dt 



fit) 



a{tf 



(33) 



-5^t )^a(t f + 5^(t ) d -^a(t f 



dt' 



dt' 



where the condition (31) is used. Since 5<fr is independent, (33) implies d<p^ t (t, x)/dt = 9<^(t, x)/dt. Together with 
the boundary condition </>^ (t, x) = ^>~ t (t, x) (recall that <^ and <$>~ x denote two different branches of the same path 
<I> C ;), one sees that the path corresponding to the stationary phase must obey 



(34) 



for all times in the rcgion(— 00, t), because they have the same initial value data. The remaining terms in (33) vanish 
since we consider FRW spacetimes with a(— 00) = and thus (33) is satisified. Thus, requiring that the phase to be 
stationary with respect to the variations of the boundary variable in (28) implies the equality of the + and — paths. 
On the other hand, independent variations implies the same condition for </> c /: 



6S[<j> c 



0. 



(35) 



i.e. <pd must obey the classical equations of motion. 

One may see that <I> C ; satisfies all the weak boundary conditions (30) and only the two of the strong boundary 
conditions (29). Therefore, if one insists on imposing the strong boundary conditions then d<j> c i(t,x)/dt = must also 
be satisfied. 

Till now we have not yet worked out the vacuum wave-functionals in (28) which directly affects the integrations over 
± (— 00). The vacuum wave-functional of the free theory is a Gaussian centered around <j) — 0. The exact form of the 
vacuum wave-functional in an interacting theory is not known but in perturbation theory there must arise corrections 
to the Gaussian functional. In any case, if the theory is expanded around <p = 0, the vacuum wave-functionals are 
expected to be oscillating wave functionals around that point. Consequently, the stationary phase approximation 
applied to ± (— 00) integrations implies 



4>d{-oo,x) = 0. 



(36) 



If one perturbatively expands around a different vacuum, say cfr — <j>o corresponding to a minimum of the potential, 
then (36) must be replaced by 4> c i{— 00, x) — (f>o- 
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We thus conclude that any classical configuration obeying the equations of motion (35) and the asymptotic boundary 
condition (36) gives rise to a stationary phase of the integral (28) (with an additional constraint d(f> c i(t,x)/dt = if 
strong boundary conditions are imposed). Here, there only appears a single boundary condition (36) and the situation 
is obviously different than the boundary conditions arising in a typical in-out path integral having two asymptotic 
regions. This is the point (1) mentioned in the introduction and it will be crucial in constructing classical solutions 
as stationary phases. 

Let us now expand the generating functional (5) around a solution <j) c i. We define new integration variables as 

<t> + = (t>cl+4> + , 

<t>~ = <t>d + 4>~ , (37) 
<fi = 4>b + <f>, 



where (fib is the boundary value of <fi c i 



4> b (x) = 4> c i(t,x). (38) 



In these new variables (5) becomes 

Z{J+,J-] mst = f D$ f n^+2?0-e 4S ^ + ^ + ^ + ]- 4S ^ + ^' J ^vl/ o [ ( / ) +(_oo)]vl/*[ ( / ) -(-oo)] 

J J -oo 

= cxp i f (j, cl (J+ - J~) J D4> j *]^2?0 + 2?^e^[^^ + ^ + ]-^^ i; ^' J_ ]*o[0 + (-oo)]*5[0-(-oo)], (39) 

— oo 

where ^ denotes the new vacuum wave-functionals and the hatted integration variables must obey 

<fi + (t , x) = <fi~ (t , x) = <f>(x) . (40) 
The new action S contains quadratic and higher order powers of the field variable, which can be written explicitly as 

S[$, J\ = ~\j d^x^—g [V^VV + V(j>) - 2J0] , (41) 

where the new potential V is given by 

V-(^) = V{j> + <fi cl ) - V(<fi cl ) - V>(cj> cl )4>. (42) 

The terms linear in (fy^ cancel out after an integration by parts since <fi c i obeys equations of motion and surface terms 
vanish owing to the boundary conditions (40). 

Differentiating with respect to J + and setting all external sources to zero one finds from (39) that 

<<fi>=4>cu (43) 

where we assume that <fi = is a minimum of V{<fi) and the path integral in (39) is perturbatively evaluated around this 
vacuum implying < <fi >= 0. A non-zero vacuum expectation value such as (43) cannot be generated in perturbation 
theory. 

It is very crucial to note that in expanding the action functionals in (35) around <fi c i, the same zeroth order term 
S[(fi c i] cancels each other in the exponential. As a result, for the in-in path integral there is no need to assume S[<f> c {\ to 
be finite to have a well defined expansion around (fid . This last property, which is also mentioned in the introduction, 
allows a more general set of field configurations to become stationary phases. 

As one would expect, in the following we will assume <fi c i to depend only on time, 

<fid = 0ci(*'). (44) 

which is suitable for cosmological applications. To warm up for our actual construction, we start studying simple 
models in flat space. 9 



9 In flat space the second line of (33) does not vanish identically. One should impose an extra "Dirichlet" boundary condition at infinity 
to set 5<^ ± (— oo) = 0. 
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= can be solved as <p c i = c\ + c 2 t'. The only solution 



Free massless scalar in flat space: The field equation 
obeying (36) is <j) c i = and thus there is no solution. 

Free massive scalar in flat space: The equation of motion </> c ; + m 2 (j> c i 
we first keep to finite and thus the solution becomes 4> c l = csm(m(t' - 
defined and thus there is no solution for this case either. 

Massless Xcj) 4 theory in flat space: The equation <f> c i + X(^ cl = has two oscillating solutions which can be expressed 
in terms of elliptic functions. The solution obeying the necessary boundary condition at to (once more keeping t 
finite first) is 



= has oscillating solutions. To impose (36) 
to))- However, to — > — oo limit is not well 



4>cl = ±c 



1/4 



sn 



1/4 



(t'-t ),-l 



(45) 



where sn is the Jacobi elliptic function, but to — > — oo limit is ill defined since the argument of sn does not converge. 
As another try, one may scale the constant c to rewrite the solution as 



*c = ±- ( 2 ~ 



1/4 



sn 



to V 2 



1/4 



(t'-f ),-l 



(46) 



This time, the argument of the sn function has a well defined limit as to — > — oo, but one ends up with the trivial 
solution 4> c i = due to the extra factor of to appeared in front. 

Although it is possible to consider more examples, these are enough for one to convince himself that in flat space 
no solution can be found for scalar fields. The reason is that the classical solutions are necessarily oscillating for all 
times and (36) only fixes the phase of the oscillation. Therefore, the solution automatically becomes ill defined at 
time t due to the infinite amount of oscillations it had performed on the way. It is clear that the expansion of the 
universe, which shows up as a cosmic friction term in the scalar field equation, might change this situation. Let us 
therefore consider some examples in dc Sitter space. 

Free massless scalar in de Sitter space: The equation <j> c i + 3Hcf) c i — can be solved as <pd = c\ + c 2 cxp(— 3Ht'). 
To impose (36) let us first keep to finite again, which gives (f> c i — </>o(l — exp[— 3H(t' — to)])- Sending to — > — oo, one 
finds c/) c i — > (fio for all finite t', which is a trivial solution. 

Free massive scalar in de Sitter space: The equation of motion <j) c i + 3H(f> c i + m 2 (fi c i = can be solved by assuming 
4>d = exp(/3t') where 1 + 3H/3 + m 2 = 0. The solution obeying (36) can be written as <fi c i — c(exp[/3 + (t' — to)] — 
exp[/3_(t' — to)])- The real part of (i is always negative and therefore (f) c i ~ > as t — > — oo for all finite t', which shows 
that there is no suitable solution. 

Massless Xcj) 4 theory in de Sitter space: In this case, the equation for cj> c i becomes 4> c i + 3H<fr c i + \<t>% = 0, which 
cannot be solved analytically. However, the dynamics is well understood: the Hubble term damps the oscillations 
about the equilibrium point <p — of X(j> 4 potential. Assuming initially that the expansion rate is small compared to 
the frequency of the oscillations, one can write <p = Af, where A represents the slowly varying amplitude and / is the 
rapidly oscillating function. Taking A » H, A <~ AH and / 3> Hf correspond to an expansion that is slow compared 
to the oscillations. Under these assumptions, the field equation can be approximately solved by imposing 



2A + 3HA = 0, 



f + XA 2 f = 0, 



(47) 



which shows that the amplitude decreases exponentially A = Aq exp(- 
elliptic function sn. Therefore, the solution can be written as 10 



-3H(t — to)/2 and / is given by the Jacobi 



<t>d ^ A 



1/4 



sn 



1/4 



J A(t')dt',--l 



(48) 



and it exponentially collapses to zero as t — > — oo. It is easy to see that if the expansion rate is not small compared 
to the initial oscillation frequency, the solution vanishes more rapidly in the limit. 

One can have a physical understanding of these negative results as follows. As noted before, the in-in path integral 
is given by the sum over all paths extending from t = — oo to t and then back to t again. The trivial path </> = 
corresponding to the " ground state" is clearly a stationary phase of the integral and the perturbation theory works 
around this path. Consider now the paths with 4>(to) = having arbitrary initial velocities at to (such paths exist in 



The period of the sn[x, —1] is given by the complete elliptic integral of first kind K(x) as 2K(— 1) ~ 2.62, which shows that the amplitude 
A in (48) also fixes the frequency of the oscillations. 
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FIG. 2. An example of a scalar potential supporting instanton-like solutions. 



the function space provided they satisfy the boundary conditions we discussed in the previous subsection) . These paths 
can be viewed as spontaneous quantum fluctuations around the vacuum and among them only stationary phases have 
a chance to give significant contribution to the in-in path integral. In fiat space, all classical solutions with <p(to) = 0, 
i.e. all vacuum fluctuations, necessarily oscillate indefinitely and their contribution to the path integral averages out 
zero, which explains the absence of solutions in that case. On the other hand, in dc Sitter space fluctuations starting 
with finite initial velocities at to = — 00 are ai l damped by the expansion of the universe near infinity so that they 
all vanish "before" escaping from the asymptotic region. Therefore, the problem is that the solutions either oscillates 
too much (flat space) or damped too much (de Sitter space). 11 

The above comments suggest where to look for nontrivial stationary phases solutions. Take a scalar field propagating 
in an expanding FRW space with a potential pictured in Fig. 2. A stationary phase (j) c i must obey 

(x d\f 

4>cl + 3-0d + 1T7 = 0. (49) 
a oq> 

Assume that we are expanding around the vacuum <j> a — 0. Solutions with small initial velocities are damped around 
4> a = and the corresponding solution in the limit to — > — oo is the trivial one <j> c i — <fi a — 0. Consider now the paths 
with larger initial velocities <p c l{to) > such that the scalar jumps over the first bump in the right and performs 
damped oscillations around the second vacuum fa- It is clear that in the to — > — oo limit the solution becomes 12 

<t>d = <f>b- (50) 

Similarly, if 4> c i{to) is chosen large enough so that the scalar overshoots the vacuum 4>b and starts performing damped 
oscillations around the locally stable vacuum <f> c , then the limiting solution corresponding to all such initial data 
becomes 

<i>d = 4>c- (51) 

Thus, all solutions starting from (/> c i(to) = with finite initial velocities asymptotically becomes either </> c ; = <j> a = 0, 
(pd = <Pb or (f> c i — <p c , which are the nontrivial stationary phases. Note that these solutions also satisfy the strong 
boundary conditions (29). 

The existence of these solutions implies that even one starts from the vacuum around <p a — 0, i.e. even when the 
vacuum wave-functionals in (5) are chosen accordingly, the generating functional (5) will get contributions from the 
vacua (f>b and <p c via the above instanton-like solutions and one should write 

Z = Z a + Z b + Z c , (52) 

where Z a , Z& and Z c can be calculated using perturbation theory around the corresponding vacuum. Moreover one 
finds 

< (f> >=(f> a + <i>b + 0c, (53) 

thus these solutions change the vacuum expectation value of the scalar from its naive perturbative value. This is very 
similar to tunneling to the perturbatively inaccessible vacua. 



As a side comment let us mention that we are focusing on solutions which has finite velocities d<f> c i(to)/dt' so that the second line of (33) 
is satisfied. It would be interesting to study solutions with infinite initial velocities still satisfying (33), which might beat the infinite 
damping of the cosmic expansion. 

Note that to — ► — oo limit is equivalent to t' — > +oo limit. 
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C. The case of finite to 

For finite to, there appears additional technical difficulties. One of the main problems is that in the absence of an 
asymptotic region even the free vacuum cannot be uniquely determined in an expanding universe [55] (recall how one 
defines the Bunch-Davies vacuum). On the other hand, in an interacting theory one cannot use the trick of giving 
a small imaginary piece to time coordinate to project onto the exact vacuum and thus even perturbation theory 
might become difficult to apply. Moreover, as we will see in a moment, the existence of stationary phases depends on 
whether one imposes the strong or the weak boundary conditions, discussed above. 

We can bypass some of these difficulties since we are employing a semiclassical approximation. Taking = 
as the minimum of the scalar potential, it is reasonable to assume that the vacuum wave-functionals are oscillating 
functionals of the field variables around this minimum. Therefore, in the stationary phase approximation, the presence 
of the vacuum wave-functionals in the path integral implies </>i(to) — as a result of ^(to) integrals in (28). The 
variation of the phase in the path integral is still given by (33) (with finite to) and setting the coefficient of the first 
line to zero again implies </>^ = Therefore, stationary phases must obey 

4>ci = 4>cl = <!><*, 4>d{to)=0, (54) 

similar to t = — oo case. 

The rest of the discussion depends on which boundary conditions arc imposed in the function space. If one assumes 
the weak boundary conditions (30), then the only way to set the second line of (33) to zero is to impose d<j> c i{to) / dt' = 
(note that 5^ (to) variations are independent), which then gives 4> c i = 0. Thus, there is no solution for weak boundary 
conditions. 

On the other hand, to satisfy the second strong condition involving a(to) in (29) one should impose Dirichlet or 
Neumann conditions for + and — branches. When the Neumann condition is chosen either for <f> + or <f>~~ , then one 
gets (j) c i = 0. Therefore, instanton-like solutions exist only when Dirichlet conditions are imposed, i.e. </' ± (to) = 0. In 
that case the second line of (33) is also satisfied since 84^ (to) = and as a result one finds 

Mh)=0, ^7^=0, (55) 

where the second condition follows from the first set of strong boundary conditions in (29). 

It is not difficult to construct nontrivial classical solutions 13 (actually infinitely many of them) satisfying (55). As an 
example consider again massless A</> 4 theory in flat space. The solution (45) already satisfies the condition 4> c i(to) = 
and imposing the second one in (55) gives 

^rt^n^.n r( 2n + i)jf(_i) 



,S/7 



t-t 



-(t'-t ),-l 



(56) 



where K(x) is the elliptic integral of first kind, K{—1) equals the half period of the Jacobi elliptic function sn[x, —1] 
and n is an integer. Note that as n getting larger, both the amplitude and the frequency of the oscillations grow, 
which shows the necessity of imposing an upper limit (cutoff) for n. The vacuum expectation value of the scalar still 
vanishes due to <j> — > —(f) symmetry (note ± signs in (56)): 

0. (57) 

Note also that as to — > — oo, (f) c i — > consistent with our earlier considerations. 
Expanding the theory around one of these solutions gives the potential 

V = 3A<^V) 2 + 2A</#V) ^ + ^ 4 , (58) 

where <j> is the fluctuation field. The potential contains time dependent coupling constants and specifically a time 
dependent mass term. Since (jJ^f 1 oscillates about zero, the shape of the potential also changes in time due to sign 
flips of the cubic term. From (56), the mass term becomes independent of A and the cubic interaction term has the 
strength \/A, which cannot arise in any perturbative expansion in A. 

Similar solutions can be seen to exist for different scalar potentials and for scalars propagating in an expanding FRW 
universe, although it might not be possible to write down analytical expressions. For instance, in the massless A^ 4 
theory defined in an exponentially expanding spacetime (which is not de Sitter space since to is finite), the solution 
become very much like (48), where the initial value of the time dependent amplitude must be quantized to satisfy 
(55). Indeed, (48) becomes more and more reliable for larger amplitudes and the solution approaches to (56). 



13 Note that the function space with strong boundary conditions is smaller than the function space with weak boundary conditions. 
Therefore, one should not surprise to see that the same path integral has stationary phases in the first space but not in the second one. 
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IV. CONCLUSIONS 



Path integral formulation gives valuable non-perturbative information about the structure of quantum field theories, 
which is otherwise hard, if not impossible, to acquire. Likewise, it is not going to be surprising to see that in-in path 
integrals will reveal some non-perturbative aspects of quantum contributions to cosmological correlations. Despite 
its evident importance, the path integral formalism applied to in-in correlation functions is not studied too much. 
Although one may encounter earlier works like [3, 5], (to our knowledge) even the very basic result on the equivalence 
of the operator and the path integral approaches in perturbation theory has been proved relatively recently in [1]. 
Similarly, a quantitative in-in path integral treatment of the cosmological perturbation theory has been given about 
two years ago [60], which is very novel. Our work confirms that in-in path integral formalism will be an important 
technique in searching for non-perturbative effects for cosmological correlation functions and some surprising results 
may still emerge. 

Including the contributions of non-trivial stationary phases to in-in path integrals, some perturbative quantum 
results are destined to change and this might have important theoretical and observational consequences. For example, 
in discussing the effects of symmetry breaking and Higgs mechanism in cosmology, the existence of instanton-like states 
must definitely be considered since they allow all locally stable minima to contribute the cosmological correlations. 
Namely, even a gauge symmetry is broken by a non-zero vacuum expectation of a scalar in a certain cosmological era, 
correlation functions might still be modified by the presence of the symmetric vacuum via instanton-like states. On 
the other hand, if time dependent solutions (i.e. instanton-like solutions whose time dependence survives the cosmic 
damping) exists similar to (56), they might affect the spectrum of cosmological perturbations (like by changing the 
tilt of the spectrum) since the spectrum is sensitive to shape of the scalar potential and these solutions give rise to 
potentials with time dependent coupling constants like (58). 

It is possible to extend this work in different directions. It would be very interesting to apply the present arguments 
to a quantum mechanical problem to test the validity of the stationary phase approximation in calculating in-in path 
integrals. As an inevitable extension, one can consider gauge fields coupled to scalars and work out how the existence 
of instanton-like solutions might have an impact on the Higgs mechanism on cosmological scales. It would also be 
an important exercise to include the gravitational degrees of freedom and determine the semiclassical gravitational 
contributions, which might affect the standard cosmological perturbation theory. Of course the problems involving 
the gauge and the gravitational fields have an extra complication related to gauge degrees of freedom, but it is worth 
to understand how the standard picture will be modified. In that way one might be more confident about inflation 
or find a caveat in its main arguments. 
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